Non-ohmic critical fluctuation conductivity of layered superconductors in magnetic 

field 
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Thermal fluctuation conductivity for a layered superconductor in perpendicular magnetic field is 
treated in the frame of the self-consistent Hartree approximation for an arbitrarily strong in-plane 
electric field. The simultaneous application of the two fields results in a slightly stronger suppression 
of the superconducting fluctuations, compared to the case when the fields are applied individually. 
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The high-temperature superconductors (HTSC) show 
a more pronounced effect of fluctuations in the normal- 
superconducting transition region, due to their high crit- 
ical temperature, small coherence length, and quasi-two- 
dimensional nature. Outside the critical region above 
T c , in the absence of magnetic field and for small elec- 
tric fields, the fluctuation conductivity can be explained 
by the Aslamazov-Larkini theory, subsequently extended 
by Lawrence and Doniach (LD)S for two-dimensional lay- 
ered superconductors, a situation very much resembling 
the crystal structure in cuprates. 

The fluctuation conductivity may be however calcu- 
lated in the linear-response approximation only for suffi- 
ciently weak electric fields, that do not perturb the fluc- 
tuation spectrum^ Reasonably high values of the elec- 
tric field accelerate the fluctuating paired carriers to the 
depairing current, and suppress the fluctuation lifetime, 
leading to deviation from the Ohm's law. In connection 
with the low-temperature superconductors, the nonlin- 
earity has been studied theoretically for the isotropic 
case^*^ while more recently Varlamov and Reggian: 6 
and Mishonov et al£ addressed the issue of the non- 
ohmic fluctuation conductivity for a layered supercon- 
ductor in an arbitrary electric field. The above men- 
tioned theories describe however the fluctuations as non- 
interacting Gaussian ones, neglecting the quartic term in 
the Ginzburg-Landau (GL) free energy. This approxi- 
mation holds not too close to the transition point, but it 
breaks down in the critical region, for higher densities of 
fluctuation Cooper pairs. 

The Ohmic fluctuation conductivity in the presence of 
a magnetic field was also initially treated^ in the Gaus- 
sian fluctuation approach, which predicted a divergence 
at T C (H) that is however not observed. The physical rea- 
son is the motion of vortices providing dissipation and 
hence a finite flux-flow conductivity. Ikeda et alw& and 
Ullah and Dorsey (UD)ii showed that the theoretical di- 
vergence can be eliminated by using the Hartree approx- 
imation, which treats self-consistently the quartic term 
in the GL free-energy expansion. This approach was ap- 
plied for the longitudinaliMi and Hall conductivity^ for 
the Nernst effect and the thermopowerjii in the linear- 
response approximation, for a layered superconductor 
under magnetic field. 

In this paper we shall treat, in the self-consistent 



Hartree approximation, the thermal fluctuation conduc- 
tivity for a layered superconductor in a perpendicular 
magnetic field and for an arbitrarily strong in-plane elec- 
tric field, a topic that, to our present knowledge, has not 
been treated yet. While the effect of the interacting fluc- 
tuations under applied magnetic field was investigated in 
the Hartree modeiifi*ii only in the linear response approx- 
imation, that is for infinitesimally small electric fields, the 
non-linear conductivity for a layered system under arbi- 
trarily strong electric field was derived in the Gaussian^ 
as well as in the Hartree approximation^ but only in the 
absence of magnetic field. 

For our purpose, we shall adopt the Langevin ap- 
proach to the time-dependent Ginzburg-Landau (TDGL) 
equation4iii The star ting point of the formalism is 
the same as in Refs. If fl and Il2l and consists of the 
LD expression of the GL free energy for a system of 
Josephson coupled superconducting planes, the gauge- 
invariant rclaxational TDGL equation governing the crit- 
ical dynamics of the superconducting order parameter, 
and the Langevin white-noise forces that model the 
thermodynamical fluctuations, satisfying the fluctuation- 
dissipation theorem. Keeping the same notations as in 
Ref. Il2l the TDGL equation for the superconducting 
order parameter in the l-ih plane will write 
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Here m and m c are effective Cooper pair masses in the 
a6-plane and along the c-axis, respectively, s is the dis- 
tance between superconducting planes, and the pair elec- 
tric charge is e§ = 2e. The order parameter has the 
same physical dimension as in the 3D case, and SI units 
are used. The perpendicular magnetic field B is gen- 
erated by the vector potential in the Landau gauge, 
A = (Q,xB,Q), with x and y the in-plane coordinates, 
and the magnetization is neglected. The GL poten- 
tial a — age is parameterized by ao = ft 2 /2m£g = 
h 2 /2m c $ c and e = ln(T/T ), with T the mean-field 
transition temperature, while £0 and £oc are, respec- 
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tively, the in-plane and out-of-plane coherence lengths 
extrapolated at T = 0. The relaxation rate T in the 
TDGL equation is given 

The Langevin white-noise forces (i (x, t) are correlated 
through (G(x,t)Cj!(x / ,f)) = 2r - 1 fc B T(5(x - x')<5(* - 
/ s j where <$(x — x') is the 2-dimensional delta- 
function concerning the in-plane coordinates. The elec- 
tric field E is assumed along the x-axis, generated by 
the scalar potential tp — —Ex. In the chosen gauge, the 
current density operator along the x direction in the 2-th 
plane will give, after averaging with respect to the noise 



as (Cq{n,k,w)q, (n',k',oj')) = 2T^ 1 k B T{2n) 3 S(k - 
k')S(q - q')S(u> - u/)8 nn >. 

By solving Eq. (0 , and taking into account the expan- 
sion form Q , one obtains further the correlation function 
of the order parameter: 

(tin (x, y, t) (x', y, <)) = 2r - 1 fc B T (7) 
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As mentioned, the quartic term in the thermo- 
dynamical potential will be treated in the Hartree 
approximation^^ by replacing the cubic term b \ipi\ 2 ipi 

in Eq. |Q with b (^ip^ 2 ^ ipi. This results in a linear 
problem with a modified (renormalized) GL potential 
a = a 
ture 



■u n ,[x' + —)(M*-M)2, 



where the expression (M* • M)~ n , is to be understood as 
the element of the inverted matrix. 

It is hereafter more convenient to rescale the integra- 
tion variables to the new ones 
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b {\^\ ), and a renormalized reduced tempera- and tQ introducc the reduced fieM magnitudes 
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We shall further introduce the Fourier transform with 
respect to the in-plane coordinate y, the layer index I, 
and time t, respectively, and also the Landau level (LL) 
representation with respect to the ^-dependence, through 
the relation: 
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where E = 16 \fik B T / 7re£o is defined as in Refs 
and 0- One can thus use further instead of the matrix 
M the following matrix M' = a ' 1 M having no physical 
dimension, so that 



ipi(x,y,t) = / — 
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where the functions u n (x) with n € N build the or- 
thonormal eigenfunction system of the harmonic oscil- 
lator hamiltonian, so that (—h 2 d 2 + e 2 l B 2 x 2 ) u n (x) = 
heoB (2n + 1) u n (x). 

The equation Q can be put further, after applying the 
expansion fll|. in the matrix form: 
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where the symmetrical tridiagonal matrix M has the non- 
zero elements: 

i f Ek \ „ he B h 2 (l — cosqs) 
' u ) + a + 
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M n+1 , n = M„, n+1 = -iT^E^ ^ Vn + 1 , 



and where the new noise terms £ g (n, k, ui), correspond- 
ing to the expansion rule l@J), are delta-correlated such 



After writing Eq. J7J) in the newly introduced variables 
and performing the integral over k' , we are able to com- 
pute the fluctuation Cooper pair density (\ip\ 2 ^ and write 

the self-consistent equation © for the renormalized re- 
duced temperature parameter e under the form: 

e=inL +gTAh r m r g-E^M > 

where we have introduced the hcrmitian ma- 
trix Q = (M'* • M') _1 and the parameter g = 
2/xoK 2 e 2 Co^s/ (nh 2 s). The expression of the quar- 
tic term coefficient 11 b = iiQH 2 e^h 2 /2m 2 was also 
taken into account, with k being the Ginzburg-Landau 
parameter k = Ao/£o- It must be noticed that the 
equation for the renormalized reduced temperature, e, is 
highly nonlinear, since the parameter e enters the M 00 
expression (|1(JH and therefore the Q matrix elements. 



3 



Analogously, starting from the correlation function Q 
written in the new variables from Eqs. © and we 
can eventually find a simpler form for the fluctuation 
conductivity: 



2/i 3 / 2 
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1 Im (Q„+i,„) 
(12) 



where the recursive properties of the «„ (x)functions and 
the Q matrix hermiticity were also used. 

In both expressions Eqs. and 112fl . the sum 

over the Landau levels must be cut-off at some maxi- 
mum index N c , reflecting the inherent UV divergence 
of the Ginzburg-Landau theory. The classical^i^ pro- 
cedure is to suppress the short wavelength fluctuating 
modes through a momentum (or, equivalently, kinetic 
energy) cut-off condition, which, in terms of the Lan- 
dau level representation write g n i 14 (he^Bjm) (n + i) < 
ca = cft 2 /2m£g, where c is an adimensional cut-off pa- 
rameter of the order of unity. A total energy cut-off 
was also proposed, 15 which eliminates the most energetic 
fluctuations and not only those with short wavelengths, 
and whose physical meaning was recently shown to follow 
from the uncertainty principle^ which imposes a limit 
to the confinement of the superconducting wave function. 
However, in the critical fluctuation region the two cut-off 
prescriptions almost coincide quantitatively, due to the 
low reduced-temperature e with respect to c, so that we 
shall apply for simplicity the cut-off procedure in its clas- 
sical form. In terms of the reduced magnetic field h, it 
writes thus h (N c + |) = c/2. In this way, the matrices 
M, M' and Q are truncated at N c + 1 lines and columns, 
and the formulae i|llfl and (|12|l can be computed numer- 
ically. 

The renormalization procedure consists thus in using 
the reduced temperature parameter e, renormalized by 
solving Eq. in the conductivity expression (|12|l . 

This procedure causes the critical temperature to shift 
downwards with respect to the bare mean-field transi- 
tion temperature To. In analogy with the Gaussian fluc- 
tuation case, we shall adopt as definition for the criti- 
cal temperature T C (E,B) the vanishing of the reduced 
temperature, e = 0. In practice, one knows experimen- 
tally the actual critical temperature T c (0,0) = T c q mea- 
sured at very low electrical field and with zero magnetic 
field, so that one would have to consider Eq. in 
the zero-fields limit. The relationship between T c q and 
To has been however already found in Ref. 12' and writes 

2gT c0 



To = T c o y/c/r + y/1 + (c/r) " . Now, having T one 

can use Eq. (|llfl for any temperature T and fields E and 
B in order to find the actual renormalized e(T, E, B), 
and further the conductivity a(T,E,B) from Eq. ltT2|) . 

In order to illustrate the main features of our model, we 
take as example a common material, like the optimally 
doped YBa 2 Cu306+£i;. Typical characteristic parameters 
are then: s = 1.17 nm for the interlayer distance, £0 = 1-2 
nm and £oc = 0.14 nm for the zero-temperature in-plane 



and out-of-plane coherence lengths, respectively, k = 70 
for the Ginzburg-Landau parameter and T c o = 92 K for 
the critical temperature under very small electric field 
and no magnetic field. We also choose for convenience 
a linear temperature extrapolation for the normal state 
resistivity which vanishes at T = 0, and has a typical 
value p N = 84 (iJJcm at T = 200 K. 

In Fig. {T^l the resistivity curves computed according 
to Eqs. (|11|) and l|12l) are shown for different magnitudes 
of the magnetic field, at a fixed value of a strong electric 
field. The zero magnetic field curve is however computed 
within the model presented in Ref. Il2l that treats the 
non-ohmic fluctuation conductivity in arbitrarily strong 
electric field and zero magnetic field, in the Hartree ap- 
proximation and with consideration of the UV cut-off. 
We can notice that the numerical results obtained for fi- 
nite magnetic fields tend to approach the curve computed 
according to Ref. while B decreases towards zero. 

Figure \T]p shows the complementary case of a fixed 
magnetic field and different magnitudes of the in- 
plane electric field. The fact that the curves ap- 
proach the result from the UD model in the vanish- 
ing electric field limit can be easily proven also ana- 
lytically. Equation (|llfl can be directly written in the 
E = limit, and since in this case the M-matrix be- 
comes diagonal, the inverted matrix Q = (M'*M') 1 
can be trivially computed as having the elements 

Q { nl = {c/ 2 + [e + r (1 - cos q') /2 + h (2n + l)] 2 } \ 
that the two integrals over ui' and q' can be immediately 
solved and yield e| B=0 = 

D„ = [{e + h + 2nh) (e + h + 2nh + r)]~ 1/2 . This is co- 
incident with the analogous equation found by UD for the 
renormalized reduced temperature in the linear response 
approximation. 

Equation l|12|) can be also processed in the vanish- 
ing electric field limit by expanding the Q-matrix up to 
the linear term in E', so that one obtains ImQ^ 1 ^ n = 

4\/6E'VhyJn+ 1 QlnQi°+ lin+1 , and after solving the in- 
tegrals over uj' and q' , the result becomes o~\ E ^ — 

(e 2 /4hs) E^=o (« + ( D n - 2^„+i + ■ This 

latter expression matches thus the linear response ap- 
proximation for the in-plane conductivity found by UD, 
as well as the corresponding formula in the paper by 
Ikeda et alM 

Figure ^ shows a further illustration of the effect ob- 
tained by the simultaneous application of the magnetic 
and electric fields, as compared to the individual effects 
produced by each field alone. The zero-fields curve is 
computed according to the LD model with the UV cut- 
off included ) 12 ! 15 namely (l6fis/e 2 ) a(e)\ E B=a = 

[e(e + r)Y 1/2 ~ [(e + c)(e + c + r 
£ + r/2)[(c + £ + r)(c + ?)r 3/2 , 
renormalized reduced temperature parame- 
ter e is given by the self-consistent Hartree 
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FIG. 1: Resistivity as a function of temperature for different combinations of the perpendicularly applied magnetic field and 
the in-plane electric field. The material parameters are given in the text. The arrows indicates the increasing direction of the 
relevant variable field, (a) Fixed electric field and variable magnetic field; the zero magnetic field curve is computed according 
to the model from Ref. Il2l (b) Variable electric field at constant magnetic field; the zero electric field curve (linear response 
limit) is computed according to the UD model. The UV cut-off parameter c = 1 in our model corresponds with limiting the 
sum on the Landau levels at the index l/2h in the UD model, (c) Effects produced by the individually (dash and dot curves) 
and the simultaneously (solid curve) applied strong electric and magnetic fields; the resistivity curve in the absence of applied 
fields (dash-dot curve) corresponds to the LD model taken in the Hartree approximation and with considering the UV cut-off. 



2gT\n (Ve + c + Ve + c + r) / L/e + 

The values for the electric and magnetic fields in Fig. 
^ were chosen such that their individual effect be almost 
similar. One can notice that the simultaneous applica- 
tion of the two fields brings however only a slight en- 
hancement of the superconducting fluctuation suppres- 
sion. Nevertheless, the non-ohmic behaviour of the fluc- 
tuation conductivity at high electric fields still remains 
quantitatively important for commonly used HTSC also 
in the presence of strong magnetic fields, so that exper- 
imental investigations could be able to discern between 
the applicability of this model in competition with the 
linear response approximation. 

In summary, we have treated the critical fluctuation 
conductivity for a layered superconductor in perpendic- 



ular magnetic field, in the frame of the self-consistent 
Hartree approximation, for an arbitrary in-plane electric 
field value. The main results are the formulae l|12|) for the 
fluctuation conductivity, and for the renormalized 
reduced-temperature parameter. In the two limit cases, 
namely for (E = 0, B > 0) and for (E > 0, B = 0), 
the corresponding solutions were found, analytically and 
respectively numerically, to reduce to previous results of 
existing theories. Qualitatively, the simultaneous appli- 
cation of the two fields results in a slight additional sup- 
pression of the superconducting fluctuation, compared to 
the case when the fields are individually applied. 
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